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Abstract 



The suppression of the individual Y(nS) states in PbPb collisions with respect to their 
yields in pp data has been measured. The PbPb and pp datasets used in the anal- 
ysis correspond to integrated luminosities of 150 ^b _1 and 230 nb -1 , respectively 
collected in 2011 by the CMS experiment at the LHC, at a center-of-mass energy 
per nucleon pair of 2.76 TeV. The Y(nS) yields are measured from the dimuon in- 
variant mass spectra. The suppression of the Y(nS) yields in PbPb relative to the 
yields in pp scaled by the number of nucleon-nucleon collisions, Raa, is measured 
as a function of the collision centrality. Integrated over centrality the Raa values are 
0.56 ± 0.08 (stat.) ± 0.07 (syst.), 0.12 ± 0.04 (stat.) ± 0.02 ( syst.), and lower than 0.10 (at 
95% confidence level), for the Y(1S), Y(2S), and Y(3S) states, respectively. The results 
demonstrate the sequential suppression of the Y(nS) states in PbPb collisions at LHC 
energies. 
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Suppression of heavy quarkonium states has been proposed as a probe of the properties of 
the hot and dense medium created in high-energy heavy-ion collisions [1J. If a deconfined 
state, often referred to as the quark-gluon plasma (QGP), is formed, the confining potential 
of heavy quark-antiquark pairs is expected to be screened because of interactions with quarks 
and gluons in the medium. The resulting dissociation of the quarkonium states depends on 
the temperature of the medium, and is expected to occur sequentially, reflecting the increasing 
values of their binding energies [2j. The Y(1S) is the most tightly bound quarkonium state, and 
is hence expected to be the one with the highest dissociation temperature. 

The prediction of the suppression pattern is complicated by various factors. These include feed- 
down contributions from higher-mass resonances into the observed quarkonium yields, as well 
as several competing nuclear and medium effects. These factors have played an important 
role in the interpretation of the charmonium measurements 0. The bottomonium family is 
expected to provide additional and theoretically cleaner probes of the deconfined medium. 
The three Y(nS) states, characterized by similar decay kinematics but distinct binding energies, 
further enable the measurement of relative state suppression, where common experimental and 
theoretical factors, and respective uncertainties, cancel. 

Measurements of the absolute Y(1S) suppression [4j and of the relative suppression of Y(2S) + 
Y(3S) with respect to Y(1S) 00 were recently reported. These analyses used PbPb (pp) data 
corresponding to an integrated luminosity of 7.3 jib~ (230 nb ) collected in 2010 (2011) at the 
same center-of-mass energy per nucleon pair of i/snn = 2.76 TeV, with the Compact Muon 
Solenoid (CMS) detector at the Large Hadron Collider (LHC). The total Y yields in PbPb (pp) 
collisions are denoted by Y|pt,pb (Y| pp ). Selecting reconstructed muons with pseudorapidity 
\t] | < 2.4 and transverse momentum pj > 4GeV/c, the Y(1S) nuclear modification factor 
Raa (defined in Eq. ^) was measured to be 0.63 ± 0.11 (stat.) ± 0.10 (syst). The double ra- 
tio Y Y(2s+3s(/Y(is) I | bFb was measure d in the same muon kinematic region to be 0.31 ^ (stat.) ± 
0.03 (syst.), indicating that the excited Y(nS) states are suppressed with respect to the Y(1S), at 
a significance of 2.4 standard deviations (c). In this Letter, an update of these measurements is 
reported, utilizing a PbPb data sample corresponding to an integrated luminosity of 150 fib -1 , 
collected in 2011 by CMS, at ,/snn = 2.76 TeV as in the previous study. This larger PbPb dataset 
together with the excellent momentum resolution of the CMS detector enables the separation of 
all three Y states below open bottom threshold for the first time in the heavy-ion environment 
and the measurement of the centrality dependence of their yields. 

A detailed description of the CMS detector can be found elsewhere [6j. Its central feature is a 
superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within 
the field volume are the silicon pixel and strip tracker, the crystal electromagnetic calorime- 
ter, and the brass/ scintillator hadron calorimeter. The silicon pixel and strip tracker measures 
charged-particle trajectories in the range \rj\ < 2.5. The tracker consists of 66 M pixel and 
10 M strip sensor elements. Muons are detected in the range \t]\ < 2.4, with detection planes 
based on three technologies: drift tubes, cathode strip chambers, and resistive plate chambers. 
Because of the strong magnetic field and the fine granularity of the tracker, the muon pj mea- 
surement based on information from the tracker alone has a resolution between 1 and 2% for a 
typical muon in this analysis. 

The CMS apparatus also has extensive forward calorimetry including two steel/ quartz-fiber 
Cerenkov hadron forward calorimeters (HF), which cover the range 2.9 < \rj\ < 5.2. These de- 
tectors are used for event selection and centrality determination in PbPb collisions. The event 
centrality observable corresponds to the fraction of the total inelastic cross section, starting at 
0% for the most central collisions and evaluated as percentiles of the distribution of the energy 
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deposited in the HF EUHl. The centrality classes used in this analysis are 50-100%, 40-50%, 
30-40%, 20-30%, 10-20%, 5-10%, and 0-5%, ordered from the lowest to the highest HF energy 
deposit. Using a Glauber-model calculation as described in Ref. (7|, the average number of nu- 
cleons participating in the collisions (Np art ) and the average nuclear overlap function (Taa) has 
been estimated for each centrality class. The Taa factor is equal to the number of elementary 
nucleon-nucleon (NN) binary collisions divided by the elementary NN cross section and can 
be interpreted as the NN-equivalent integrated luminosity per heavy-ion collision, at a given 
event centrality Q. 

The Y states are identified through their dimuon decay. The events are selected online with a 
hardware-based trigger requiring two muon candidates in the muon detectors. More stringent 
muon quality requirements are imposed in the PbPb case relative to the pp online selection. 
No explicit momentum or rapidity thresholds are applied at trigger level. For the PbPb data, 
events are preselected offline if they contain a reconstructed primary vertex comprising at least 
two tracks, and the presence of energy deposits larger than 3 GeV in at least three towers in each 
of the two HF calorimeters. These criteria reduce contributions from single-beam interactions, 
ultra-peripheral electromagnetic interactions, and cosmic-ray muons. 

Muons are reconstructed by matching tracks in the muon detectors and silicon tracker. The 
same offline reconstruction algorithm and selection criteria are applied to the PbPb and pp 
data samples. The muon candidates are required to have a transverse (longitudinal) distance 
of closest approach to the event vertex smaller than 3 (15) cm. Muons are only kept if the 
part of their trajectory in the tracker has 11 or more hits and the x 2 P er degree of freedom of 
the combined and tracker-only fits are lower than 20 and 4, respectively. Pairs of oppositely 
charged muons are considered dimuon candidates if the x 2 fit probability of the tracks origi- 
nating from a common vertex exceeds 5%. This removes background arising primarily from the 
displaced, semileptonic decays of charm and bottom hadrons. Only muons with pj > 4GeV/c 
and \rj\ < 2.4 are considered, as in Ref. 0. The dimuon pj distribution of the selected candi- 
dates extends down to zero and has a mean of about 6GeV/c. The resultant dimuon invariant 
mass spectra are shown in Fig. [I] for the PbPb and pp datasets. The three Y(nS) peaks are 
clearly observed in the pp case; the Y(3S) state is not visible above the dimuon continuum in 
PbPb collisions. 

Simulated Monte Carlo (MC) events are used to optimize muon selection cuts and to evaluate 
efficiencies. Signal Y(nS) events are generated using PYTHIA 6.424 [10|, with non-relativistic 
quantum chromodynamics matrix elements tuned by comparison with CDF data |TT|. Under- 
lying heavy-ion events are produced with the hydjet 1.6 ffl2| event generator. The detector 
response is simulated with GEANT4 [13J. The signal candidates are embedded in the underly- 
ing PbPb events, at the level of detector hits and with matching vertices. The resulting embed- 
ded events are then processed through the trigger emulation and the full event reconstruction 
chain. 

An extended unbinned maximum likelihood fit to the two invariant mass spectra shown in 
Fig. [T]is performed to extract the Y(nS) yields, following the method described in Refs. I[51H4|. 
The measured mass lineshape of each Y(nS) state is parameterized by a "Crystal Ball" (CB) 
function, i.e. a Gaussian resolution function with the low-side tail replaced by a power law 
describing final-state radiation. The mass differences between the states are fixed to their world 
average values l|T5| and the mass resolution is forced to scale with the resonance mass. In 
our previous measurement (5), the signal shape parameters were fixed from MC simulation, 
including the mass resolution and CB tail parameters. The current 20-fold larger PbPb dataset 
allows these constraints to be released, but the shape parameters are treated as common for 
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Figure 1: Dimuon invariant-mass distributions in PbPb (left) and pp (right) data at ^s^n = 
2.76 TeV. The same reconstruction algorithm and analysis selection are applied to both datasets, 
including a transverse momentum requirement on single muons of pj > 4GeV/c. The solid 
(signal + background) and dashed (background-only) curves show the results of the simulta- 
neous fit to the two datasets. 



both PbPb and pp datasets via a simultaneous fit. 

The background model for the pp dataset consists of a second-order polynomial, as was used 
in Ref. |5), while the larger PbPb dataset requires a more detailed background model. The 
px > 4GeV/c muon selection threshold causes a depletion of dimuon candidates in the lower 
part of the 7-14 GeV/ c 2 mass fitting range. The PbPb background model consists of an exponen- 
tial function multiplied by an error function describing the low-mass turn-on. The background 
parameters are determined from the fit. This nominal model accurately describes the mass side- 
bands in the opposite-sign muon signal sample, shown in Fig. [T] (left), as well as the alternative 
estimates of the shape of the combinatorial background obtained from like-sign muon pairs or 
via a "track-rotation" method. In the latter method [16| the azimuthal angular coordinate of 
one of the muon tracks is rotated by 180 degrees. 

The ratios of the observed yields, not corrected for differences in acceptance and efficiency, of 
the Y(2S) and Y(3S) states to the Y(1S) state, in the PbPb and pp data, are 

Y(2S)/Y(lS)| pp = 0.56 ± 0.13 (stat.) ± 0.02 (syst.) , (1) 

Y(2S)/Y(lS)| PbPb = 0.12 ± 0.03 (stat.) ± 0.02 (syst.) , 

Y(3S)/Y(lS)| pp = 0.41 ± 0.11 (stat.) ± 0.04 (syst.) , 

Y(3S)/Y(lS)| Pb Pb = 0.02 ± 0.02 (stat.) ± 0.02 (syst.) (< 0.07 at 95% confidence level) , 

where the systematic uncertainty arises from the fitting procedure, as described below. For the 
Y(3S) to Y(1S) ratio in PbPb, a 95% confidence level (CL) limit is set, based on the Feldman- 
Cousins statistical method ffl"7|. 

The measurement of the ratio of the Y(nS)/Y(lS) ratios in PbPb and pp collisions benefits from 
an almost complete cancellation of possible acceptance or efficiency differences among the re- 
constructed resonances. The simultaneous fit to the PbPb and pp mass spectra gives the double 
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Figure 2: Centrality dependence of the double ratio (left) and of the nuclear modification fac- 
tors (right) for the Y(1S) and Y(2S) states. The relative uncertainties from N par t-independent 
quantities (pp yields and, for the Raa, also integrated luminosity) are represented by the boxes 
at unity and are not included in the data points as these uncertainties do not affect the point- 
to-point trend. The event centrality bins used are indicated by percentage intervals. The results 
are available in tabulated form in App.|A| 



ratios 

Y(2S)/Y(lS)| PbPb 
Y(2S)/Y(lS)| pp 

Y(3S)/Y(lS)| PbPb 
Y(3S)/Y(lS)|p P 



0.21 ± 0.07 (stat.) ± 0.02 (syst.) , (2) 
0.06 ± 0.06 (stat.) ± 0.06 (syst.) (< 0.17 at 95% CL) . 



The systematic uncertainties from the fitting procedure are evaluated by varying the fit func- 
tion as follows: fixing the CB tail and resolution parameters to MC expectations, allowing 
for differences in these parameters between PbPb and pp, and constraining the background 
parameters with the like-sign and track-rotated spectra. An additional systematic uncertainty 
(1%), estimated from MC simulation, is included to account for possible imperfect cancellations 
of acceptance and efficiency. 

The double ratios, defined in Eq. are expected to be compatible with unity in the absence of 
suppression of the excited states relative to the Y(1S) state. The measured values are, instead, 
considerably smaller than unity. The significance of the observed suppression exceeds 5 a. 

In order to investigate the dependence of the suppression on the centrality of the collision, 
the double ratio Y Y(2s(/Y(iS)'| bPb * s displayed as a function of N par t in Fig. |ij (left). The results 
are constructed from the single ratio Y(2S)/Y(lS)|p b p b measured in bins of PbPb centrality, 
using the pp ratio as normalization. No pronounced centrality dependence is observed. More 
data, in particular more pp collisions, are needed to establish possible dependences on dimuon 
kinematic variables. 



Absolute suppressions of the individual Y states and their dependence on the collision central- 
ity are studied using the nuclear modification factor, Raa, defined as the yield per nucleon- 
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nucleon collision in PbPb relative to that in pp. The Raa observable, 

r aa — ^pp Y(nS)|pbPb £ pp ^ 
Taa^MB Y(nS)| pp epbPb 

is evaluated from the ratio of total Y(nS) yields in PbPb and pp collisions corrected for the 
difference in efficiencies £ pp / £pbpb/ with the average nuclear overlap function Taa, number of 
minimum-bias events sampled by the event selection Nmb, and integrated luminosity of the pp 
dataset £ pp accounting for the normalization. The centrality-integrated (0-100%) Raa values 
for the individual Y states are: 

Raa(Y(1S)) = 0.56 ± 0.08 (stat.) ± 0.07 (syst.), (4) 
R A a(Y(2S)) = 0.12 ±0.04 (stat.) ± 0.02 (syst.), 

Raa(Y(3S)) = 0.03 ± 0.04 (stat.) ± 0.01 (syst.) (< 0.10 at 95% CL) . 

As the Y(3S) peak is not observed above the dimuon continuum (significance less than one stan- 
dard deviation), an upper limit is also presented. The results for the Y(1S) and Y(2S) obtained 
by performing the measurement in ranges of centrality are displayed in Fig. [2] (right). 

Each factor entering in Eq. (j3l contributes to the Raa uncertainty, including £ pp (6%) and Taa 
(4-15%, from central to peripheral collisions). The systematic uncertainties from the fitting pro- 
cedure, used in the determination of the Y(1S) (4-9%), Y(2S) (10-40%), and Y(3S) (14%) signal 
yields, are estimated as previously described for the double-ratio measurement. The ratio of 
efficiencies in Eq. (|3j is estimated from MC simulation to deviate by less than 7% from unity 
for the centrality bins considered. Systematic uncertainties on the efficiency ratio are estimated 
by considering variations of simulated kinematic distributions (5-7%) and from differences in 
the efficiency ratio estimations from data and MC simulations (3%). For the former source, 
uncertainties are estimated by applying a weight to the generated Y pj and \y\ distributions 
that increases linearly from 0.7 to 1.3 over the ranges < pj < 20GeV/c and < \y\ < 2.4. 
For the latter source, reconstruction and trigger selection efficiencies are estimated employing 
a tag-and-probe method |4lll4"]|, using muons from ]/xp decays in PbPb and pp simulations as 
well as in collision data. 

The results indicate a significant suppression of the Y(nS) states in heavy-ion collisions com- 
pared to pp collisions at the same per-nucleon-pair energy. The data support the hypothesis of 
increased suppression of less strongly bound states: the Y(1S) is the least suppressed and the 
Y(3S) is the most suppressed of the three states. The Y(1S) and Y(2S) suppressions are observed 
to increase with collision centrality. The suppression of Y(2S) is stronger than that of Y(1S) in all 
centrality ranges, including the most peripheral bin. It should be noted that this bin (50-100%) 
is rather wide and mostly populated by more central events (closer to 50%). For this most pe- 
ripheral bin, the Y(1S) nuclear modification factor is consistent with unity, 1.01 ± 0.12 (stat.) ± 
0.22 (syst.), while for the most central bin (0-5%), R A a is 0.41 ± 0.04 (stat.) ± 0.07 (syst.) indicat- 
ing a significant suppression. The observed Y(nS) yields contain contributions from decays of 
heavier bottomonium states and, thus, the measured suppression is affected by the dissociation 
of these states. This feed-down contribution to the Y(1S) state was measured to be of the or- 
der of 50% lfT8| , however the precision of that measurement precludes quantitative conclusions 
about the suppression of the directly produced Y(1S). 

In addition to QGP formation, differences between quarkonium production yields in PbPb 
and pp collisions can also arise from cold-nuclear-matter effects ||T9|. However, such effects 
should have a small impact on the double ratios reported here. Initial-state nuclear effects are 
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expected to affect similarly each of the three Y states, thereby canceling out in the ratio. Final- 
state "nuclear absorption" becomes weaker with increasing energy [20J and is expected to be 
negligible at the LHC [21 j. Future high-statistics heavy-ion, proton-proton, and proton-nucleus 
runs at the LHC will provide further quarkonium measurements, which will help to disentan- 
gle cold-nuclear from hot-medium effects and to attain a more thorough characterization of the 
properties of the produced medium. 

In conclusion, the first observation of the sequential suppression of the Y(nS) states in heavy- 
ion collisions has been reported, in ^/s^n = 2.76 TeV PbPb collisions by the CMS experiment at 
the LHC, extending the previous CMS bottomonium measurements HH5). The Y(2S) and Y(3S) 
resonances are suppressed with respect to the Y(1S) state, with a significance exceeding 5 a. The 
nuclear modification factors for the Y(nS) states were also measured, with the individual Y(1S), 
Y(2S), and Y(3S) states suppressed by factors of about 2, 8, and larger than 10, respectively. 
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A Central ity-dependent results (tabulated format) 



A Centrality-dependent results (tabulated format) 



Table 1: Double ratio versus centrality. 



Np ar t (centrality) 


Y(2S)/Y(lS)| Pb pb i „«. t i „ „ t i o.1 n Vn1 
Y(2S)/Y(lS)| pp ± Stat ± S y St ± g loDal 


381 ( 0-5%) 
329 ( 5-10%) 
261 (10-20%) 
187 (20-30%) 
130 (30-40%) 

86 (40-50%) 
22 (50-100%) 


0.270 ± 0.146 ± 0.031 ± 0.032 
0.102 ± 0.189 ± 0.042 ± 0.012 
0.138 ± 0.108 ± 0.041 ± 0.017 
0.428 ± 0.124 ± 0.039 ± 0.051 
0.344 ± 0.138 ± 0.027 ± 0.041 
0.420 ± 0.220 ± 0.048 ± 0.050 
0.304 ± 0.154 ± 0.040 ± 0.037 



Table 2: Y(1S) and Y(2S) nuclear modification factors, Raa, versus centrality. 



Np ar t (centrality) 


Y(1S) R AA ± stat. ± syst. ± global 


Y(2S) R A a ± stat. ± syst. ± global 


381 ( 0-5%) 
329 ( 5-10%) 
261 (10-20%) 
187 (20-30%) 
130 (30^0%) 

86 (40-50%) 
22 (50-100%) 


0.411 ± 0.043 ± 0.048 ± 0.056 
0.432 ± 0.048 ± 0.046 ± 0.059 
0.484 ± 0.040 ± 0.049 ± 0.066 
0.614 ± 0.053 ± 0.075 ± 0.084 
0.681 ± 0.069 ± 0.085 ± 0.093 
0.590 ± 0.096 ± 0.086 ± 0.080 
1.005 ± 0.121 ± 0.176 ± 0.137 


0.111 ± 0.061 ± 0.020 ± 0.024 
0.044 ± 0.060 ± 0.019 ± 0.009 
0.068 ± 0.053 ± 0.022 ± 0.014 
0.260 ± 0.079 ± 0.043 ± 0.055 
0.237 ± 0.098 ± 0.039 ± 0.050 
0.251 ± 0.138 ± 0.050 ± 0.054 
0.300 ± 0.157 ± 0.069 ± 0.064 
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